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Abstract 

We introduce a new variant of the geometric Steiner arborcsccncc problem, motivated by 
the layout of flow maps. Flow maps show the movement of objects between places. They 
, reduce visual clutter by bundling lines smoothly and avoiding self-intersections. To capture 

these properties, our angle-restricted Steiner arborescences, or flux trees, connect several targets 
', to a source with a tree of minimal length whose arcs obey a certain restriction on the angle they 

■ form with the source. 

We study the properties of optimal flux trees and show that they are planar and consist of 
logarithmic spirals and straight lines. Flux trees have the shallow-light property. We show that 
computing optimal flux trees is NP-hard. Hence we consider a variant of flux trees which uses 
only logarithmic spirals. Spiral trees approximate flux trees within a factor depending on the 
angle restriction. Computing optimal spiral trees remains NP-hard, but we present an efficient 
2-approximation, which can be extended to avoid "positive monotone" obstacles. 

1 Introduction 

. Flow maps are a method used by cartographers to visualize the movement of objects between 

^sD ! places [HI [19]. One or more sources are connected to several targets by arcs whose thickness 

corresponds to the amount of flow between a source and a target. Good flow maps share some 
(f-) ■ common properties. They reduce visual clutter by merging (bundling) lines as smoothly and 
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frequently as possible. Furthermore, they strive to avoid crossings between lines. Flow trees, that 
is, single-source flows, are drawn entirely without crossings. Flow maps that depict trade often 
route edges along actual shipping routes. In addition, flow maps try to avoid covering important 
map features with flows to aid recognizability. Most flow maps are still drawn by hand and none 
of the existing algorithms (that use edge bundling), can guarantee to produce crossing- free flows. 

In this paper we introduce a new variant of geometric minimal Steiner arborescences, which 
captures the essential structure of flow trees and serves as a "skeleton" upon which to build high- 
quality flow trees. Our input consists of a point r, the root (source), and n points t\,...,t n , 
the terminals (targets). Visually appealing flow trees merge quickly, but smoothly. A geometric 
minimal Steiner arborescence on our input would result in the shortest possible tree, which naturally 
merges quickly. A Steiner arborescence for a given root and a set of terminals is a rooted directed 
Steiner tree, which contains all terminals and where all edges are directed away from the root. 
Without additional restrictions on the edge directions (as in the rectilinear case or in the variant 
proposed in this paper), a geometric Steiner arborescence is simply a geometric Steiner tree with 
directed edges. However, Steiner arborescences have angles of 2n/3 at every internal node and 
hence are quite far removed from the smooth appearance of hand-drawn flow maps. Our goal is 
hence to connect the terminals to the root with a Steiner tree of minimal length whose arcs obey 
a certain restriction on the angle they form with the root. 



*A preliminary version of this paper will appear at the 22nd International Symposium on Algorithms and Com- 
putation (ISAAC 2011). B. Speckmann and K. Verbeek are supported by the Netherlands Organisation for Scientific 
Research (NWO) under project no. 639.022.707. 
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Figure 1: Flow maps from our companion paper [7] based on angle-restricted Steiner arborescences: 
Migration from Colorado and whisky exports from Scotland. 

Specifically, we use a restricting angle a < tt/2 to control 
the direction of the arcs of a Steiner arborescence T. Con- 
sider a point p on an arc e from a terminal to the root (see 
Figure[2]). Let 7 be the angle between the vector from p to the 
root r and the tangent vector of e at p. We require that 7 < a 
for all points p on T. We refer to a Steiner arborescence Figure 2: The angle restriction, 
that obeys this angle restriction as angle-restricted Steiner 

arborescence, or simply flux tree. Here and in the remainder of the paper it is convenient to direct 
flux trees from the terminals to the root. Also, to simplify descriptions, we often identify the nodes 
of a flux tree T with their locations in the plane. 

In the context of flow maps it is important that flux trees can avoid obstacles, which model 
important features of the underlying geographic map. Furthermore, it is undesirable that terminals 
become internal nodes of a flux tree. We can ensure that our trees never pass directly through 
terminals by placing a small triangular obstacle just behind each terminal (as seen from the root). 
Hence our input also includes a set of m obstacles Bi,... , B m . We denote the total complexity 
(number of vertices) of all obstacles by M. In the presence of obstacles our goal is to find the 
shortest flux tree T that is planar and avoids the obstacles. 

The edges of flux trees are by definition "thin", but their topology and general structure are 
very suitable for flow trees. In a companion paper [7j we describe an algorithm that thickens and 
smoothes a given flux tree while avoiding obstacles. Figure [1] shows two examples of the maps 
computed with our algorithm, further examples and a detailed discussion of our maps can be found 

in m- 

Related work. There is a multitude of related work on both the practical and the theoretical side 
of our problem and consequently we cannot cover it all. 

One of the first systems for the automated creation of flow maps was developed by Tobler in 
the 1980s [HEO]- His system does not use edge bundling and hence the resulting maps suffer from 
visual clutter. In 2005 Phan et al. [H] presented an algorithm, based on hierarchical clustering of 
the terminals, which creates flow trees with bundled edges. This algorithm uses an iterative ad-hoc 
method to route edges and is often unable to avoid crossings. A second effect of this method is 
that flows are often routed along counterintuitive routes. The quality of the maps can be improved 
by moving the terminals, which, however, is considered to be confusing for users by cartography 
textbooks [19]. Recent papers from the information visualization community explore alternative 
ways to visualize flows, by using multi-view displays [10j . animations over time [3], or mapping 
techniques close to treemaps |21j . 

There are many variations on the classic Steiner tree problem which employ metrics that are 
related to their specific target applications. Of particular relevance to this paper is the rectilinear 
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Steiner arborescence (RSA) problem, which is defined as follows. We are given a root (usually at 
the origin) and a set of terminals t±, . . . , t n in the northeast quadrant of the plane. The goal is to 
find the shortest rooted rectilinear tree T with all edges directed away from the root, such that T 
contains all points t\, . . . , t n . For any edge of T from p = (x p , y p ) to q = (x q , y q ) it must hold that 
x p < x q and yp < y q . If we drop the condition of rectilinearity then we arrive at the Euclidean 
Steiner arborescence (ESA) problem. In both cases it is NP-hard [17} [T8] to compute a tree of 
minimum length. Rao et al. |16j give a simple 2-approximation algorithm for minimum rectilinear 
Steiner arborescences. Cordova and Lee [8] describe an efficient heuristic which works for terminals 
located anywhere in the plane. Ramnath [15] presents a more involved 2-approximation that can 
also deal with rectangular obstacles. Finally, Lu and Ruan [12] developed a PTAS for minimum 
rectilinear Steiner arborescences, which is, however, more of theoretical than of practical interest. 

Conceptually related are gradient- constrained minimum networks which are studied by Brazil et 
al. [3 E] motivated by the design of underground mines. Gradient-constrained minimum networks 
are minimum Steiner trees in three-dimensional space, in which the (absolute) gradients of all edges 
are no more than an upper bound m (so that heavy mining trucks can still drive up the ramps 
modeled by the Steiner tree). Krozel et al. |llj study algorithms for turn-constrained routing with 
thick edges in the context of air traffic control. Their paths need to avoid obstacles (bad weather 
systems) and arrive at a single target (the airport). The union of consecutive paths bears some 
similarity with flow maps, although it is not necessarily crossing-free or a tree. 

Results and organization. In Section [2] we derive properties of optimal (minimum length) flux 
trees. In particular, we show that they are planar and that the arcs of optimal flux trees consist of 
(segments of) logarithmic spirals and straight lines. Flux trees have the shallow-light property [3], 
that is, we can bound the length of an optimal flux tree in comparison with a minimum spanning 
tree on the same set of terminals and we can give an upper bound on the length of a path between 
any point in a flux tree and the root. They also naturally induce a clustering on the terminals and 
smoothly bundle lines. Unfortunately we can show that it is NP-hard (Section 14. ip to compute 
optimal flux trees. Hence, in Section [3] we introduce a variant of flux trees, so called spiral trees. 
The arcs of spiral trees consist only of logarithmic spiral segments. We prove that spiral trees 
approximate flux trees within a factor depending on the restricting angle a. Our experiments show 
that a = 7r/6 is a reasonable restricting angle, in this case the approximation factor is sec(a) ~ 1.15. 
In Section [4.11 we show that computing optimal spiral trees remains NP-hard. For a special case, 
we give an exact algorithm in Section [4.21 that runs in 0(n 3 ) time. In Section [4.31 we develop a 2- 
approximation algorithm for spiral trees that works in general and runs in O(nlogn) time. Finally, 
in Section [5] we extend our approximation algorithm (without deteriorating the approximation 
factor) to include "positive monotone" obstacles. On the way, we develop a new 2-approximation 
algorithm for rectilinear Steiner arborescences in the presence of positive monotone obstacles. Both 
algorithms run in 0((n + M) log(n + M)) time, where M is the total complexity of all obstacles. 

2 Optimal flux trees 

Recall that our input consists of a root r, terminals t±, . . . ,t n , and a restricting angle a < ir/2. 
Without loss of generality we assume that the root lies at the origin. Recall further that an optimal 
flux tree is a geometric Steiner arborescence, whose arcs are directed from the terminals to the root 
and that satisfies the angle restriction. We show that the arcs of an optimal flux tree consist of line 
segments and parts of logarithmic spirals (Property [Q) , that any node except for the root has at 
most two incoming arcs (Property [2|) , and that an optimal flux tree is planar (Property [3|). Finally, 
flux trees (and also spiral trees) have the shallow-light property (Property U|) . 
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Figure 3: Spirals and spiral regions. 



Spiral regions. For a point p in the plane, we consider 
the region 1Z P of all points that are reachable from p 
with an angle-restricted path, that is, with a path that 
satisfies the angle restriction. Clearly, the root r is al- 
ways in lZ, p . The boundaries of 1Z P consist of curves 
that follow one of the two directions that form exactly 
an angle a with the direction towards the root. Curves 
with this property are known as logarithmic spirals (see 
Figure [3]). Logarithmic spirals are self-similar; scaling a 
logarithmic spiral results in another logarithmic spiral. 

Logarithmic spirals are also self- approaching as defined by Aichholzer et al. [2], who give upper 
bounds on the lengths of (generalized) self- approaching curves. As all spirals in this paper are loga- 
rithmic, we simply refer to them as spirals. For a < tt/2 there are two spirals through a point. The 
right spiral Sp is given by the following parametric equation in polar coordinates, where p = (R, 4>): 
R(t) = Re~ l and (p(t) = <f> + tan(a)t. The parametric equation of the left spiral S~ is the same 
with a replaced by —a. Note that a right spiral S p can never cross another right spiral (the 
same holds for left spirals). The spirals S p and S~ cross infinitely often. The reachable region 1Z P 
is bounded by the parts of 5+ and S~ with < t < tt cot (a). We therefore call 1Z P the spiral region 
of p. It follows directly from the definition that for all q £ 1Z P we have that lZ q C 1Z P . 



Lemma 1 The shortest angle-restricted path between a point p and a point q 6 1Z P consists of a 
straight segment followed by a spiral segment. Either segment can have length zero. 

Proof. Consider the spirals 5+ and S~ through q, specifically the parts with t < (see Figure HJ. 
Any point on the opposite side of the spirals as p is unable to reach q. Thus any shortest path from 
p to q cannot cross either of these spirals. If we see these spirals as obstacles and ignore the angle 
restriction for now, the shortest path ir is simply a straight segment followed by a spiral segment. 
Now consider any point u on tt. Because <S+ and S^ cannot cross (same for S~ and S~), we get 
that q S 1Z U . Therefore tt also satisfies the angle restriction. □ 

Property 1 An optimal flux tree consists of straight segments and spiral segments. 

Proof. Consider an optimal flux tree T. Now replace all edges between a terminal/the root and 
a Steiner node by the shortest angle-restricted path between the two points. This can only shorten 
T. By Lemma [H the resulting flux tree consists of only straight segments and spiral segments. □ 

Property 2 Every node in an optimal flux tree T, other than the root r, has at most two incoming 
edges. 

Proof. Assume T contains a node at p with at least three incoming edges. Pick one of the 
incoming edges e that is not leftmost or rightmost and let q be the other endpoint of e. Let ei and 






Figure 4: Shortest path. Figure 5: An optimal flux tree (a = vr/6). Figure 6: Property 2. 
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en be the leftmost and rightmost incoming edges of p (see Figure E]) . Now consider the straight line 
£ from q to the root r. Assume without loss of generality that £ passes p on the left side (the right 
side is symmetric with en) or £ goes through p. We claim that we can locally improve the length 
of T by moving the endpoint at p of e along ex- The angle between e and ex at p is at most a. 
Because a < tt/2 and because locally moving the endpoint of e along ex will not make the spiral 
segment of e longer, this will shorten the tree T. Also, locally moving the endpoint of e along ex 
cannot suddenly violate the angle restriction (assuming that a > 0). Contradiction. □ 

Property 3 Every optimal flux tree is planar. 

Proof. Assume two edges e\ (from p\ to q\) and e<i (from p2 to 52) cross. Let u be the crossing 
between e\ and &i- Now simply remove the part of e\ from u to q\. There is still a connection from 
pi to r via q%, so the resulting tree is still a proper flux tree. Also, removing a segment cannot 
violate the angle restriction and makes the tree shorter. Contradiction. □ 

The last property requires a more involved proof. We postpone the proof of this property until 
Section \'6.1\ Let d T (p) be the distance between p and r in a flux tree T and let d(p) be the Euclidean 
distance between p and r. 

Property 4 The length of an optimal flux tree T is at most (3((sec(a) + csc(a)) log n) times the 
length of the minimum spanning tree on the same set of terminals. Also, for every point p £ T, 
d T (p) < sec(a)d(p). 

3 Spiral trees 

In this section we introduce spiral trees and prove that they approximate flow trees. The arcs of a 
spiral tree consist only of spiral segments of a given a (see Figure [7j). In other words, an optimal 
spiral tree is the shortest flow tree that uses only spiral segments. Spiral trees satisfy the angle 
restriction by definition. Any particular arc of a spiral tree can consist of arbitrarily many spiral 
segments. That is, any arc of the spiral tree can switch between following its right spiral and 
following its left spiral an arbitrary number of times. The length of a spiral segment can easily be 
expressed in polar coordinates. Let p = (Ri,4>i) and q = (i?2,<fe) be two points on a spiral, then 
the distance D(p, q) between p and q on the spiral is 



Consider the shortest spiral path — using only spiral segments — between a point p and a point q 
reachable from p. The reachable region for p is still its spiral region 1Z P , so necessarily q € 1Z P . The 
length of a shortest spiral path is given by Equation [TJ The shortest spiral path is not unique, in 
particular, any sequence of spiral segments from p to q is shortest, as long as we move towards the 
root. 






Figure 7: Spiral tree with spiral regions. 



Figure 8: T and T". 
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Theorem 1 The optimal spiral tree T' is a sec(a) -approximation of the optimal flux tree T. 



Proof. Let Cr be a circle of radius R with the root r as center. A lower bound for the length of 
T is given by L(T) > J °° |TD Cr\cIR, where ITPiCrJ counts the number of intersections between 
the tree T and the circle Cr. Using Equation [TJ the length of T' is L(T') = sec(a) j Q \T' n Cr\cIR. 
Now consider the spiral tree T" with the same nodes as T, but where all arcs between the nodes 
are replaced by a sequence of spiral segments (see Figured]). For a given circle Cr, this operation 
does not change the number of intersections of the tree with Cr, i.e. \T n Cr\ = \T" n Cr\. So we 
get the following: 



Next to the fact that optimal spiral trees are a good approximation of optimal flux trees, they also 
maintain important properties of optimal flux trees, namely Properties [2] and [3j 

Lemma 2 An optimal spiral tree is planar and every node, other than the root, has at most two 
incoming edges. The root has exactly one incoming edge. 

Proof. First of all, only two spirals go through a single point: the left and the right spiral. So 
every node other than the root r has at most two incoming arcs, otherwise there is a repeated spiral 
segment which can be removed. Furthermore, the same arguments as in the proof of Property [3] 
yield that also the optimal spiral tree is planar. □ 

When we approximate an optimal flux tree by a spiral tree, we can further reduce the length of 
the tree by replacing every arc by the shortest angle-restricted path between its endpoints. This 
operation does not improve the approximation factor, but it improves the tree visually. 

3.1 Shallow-Light Property 

In the following we prove the shallow-light property for optimal spiral trees. That is, we bound the 
length of an optimal spiral tree in comparison with a minimum spanning tree on the same set of 
terminals and we give an upper bound on the length of a path between any point in a spiral tree 
and the root. Since for flux trees the length of such paths and the total length are not larger, we 
can conclude that optimal flux trees also have the shallow-light property (Property . The second 
part of the property (shallowness) is easy to see. The path of a node to the root in any spiral tree 
is by Equation [1] bounded by sec(a) times the distance of the node to the root. 

We now show that the length of an optimal spiral tree approximates the length of a minimum 
spanning tree by a factor of O((sec(o) + csc(a)) log n). We build a spiral tree in the following way. 
First we find a short cycle through the points. We then take a matching based on this cycle and 
pairwise join points by spiral segments. This results in \n/2] components. On these we again find 
a matching and pairwise join them and so on. We need to ensure that the set of spiral segments 
used in the construction is compatible with a spiral tree. 

Throughout this section we will assume that the root of the tree is placed at the origin. We call 
a sequence of spiral segments between two points inward going if the distance from the segments to 
the origin have no local maximum except possibly at the two points. In particular, if a point is in 
the spiral region of another, a path of decreasing distance to the root from the outer point to the 
inner one would be inward going. Any pair of points can be joined by an inward going sequence 
of two spiral segments and inward going sequences are compatible with spiral trees if we use the 
point with smallest distance to the root as join node. 




□ 
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We need to bound the length of such a sequence. For this we first bound the length of a spiral 
segment. Let p, p' be two points on a spiral segment with polar coordinates p = (R,(j>) an d 
p' = (R',(f)'). Equation [T] gives us a bound in terms of R,R'. To bound the length in terms of 
R, <f>, (j)' , let us assume R' < R. The other case is analogous. We consider the parametric equations 
of the spiral through p and p' with (-R(O), (f)(0)) = (R, (j)). For p' we obtain the equations R' = Re~ l 
and (j)' = (j) + tan(a)t (or (j)' = (j) — tan(a)t depending on whether the points lie on a left or right 
spiral). Solving for R' yields R' = Re~^ - ^D/ tan ( a ). Inserting this into Equation Q] gives 

D(p,p) = sec{a)R(l - e -(l*'-«/W«)) . (2) 

Equation [2] has several consequences. Given two points p, q that do not lie in the spiral regions of 
each other. Assume we have two sequences of spiral segments, each connecting p and q such that 
no ray through the origin intersect a sequence twice, and such that parameterized by the angle 4> 
to the origin, one sequence has a smaller or equal distance to the origin for all (j). Then sweeping 
over the angle and summing up the contributions of Equation [2] gives that the sequence closer to 
the origin has a smaller (or equal) total arc length. Thus, the shortest connection between p and 
q is obtained by simply joining p and q by an inward going sequence of two spiral segments. 

Another consequence of Equation[2]is the following. Again p = (R, </>), q = (R 1 , (j)') are two points 
that do not lie in the spiral regions of each other. Further assume R, </>, (p' are given, but for R' we 
only know R' < R. Then the arc length of the inward going sequence of two spiral segments joining 
p and q (using the angle range between <j) and (j)') is maximized for R = R' . This follows from 
the same argument as above, i.e., the resulting sequence of spiral segments dominates all others in 
terms of distance to the origin. 

So far we have not linked the arc length of the spiral segments between two points with the 
Euclidean distance between the points. We do this by the following lemma. 

Lemma 3 Two points in the plane of distance D can be connected by an inward going path of 
logarithmic spirals of angle a such that the summed length of the spiral segments is bounded by 
3D max(sec(a), csc(a)). The path uses at most two spiral segments. 

Proof. Let p\, p2 be two points of distance D with polar coordinates p\ = (Ri,4>i) and p2 = 
(i?2, 02 )• Without loss of generality we assume that R± < R2, 4>\ < 4>2 and 4>2 — 4>\ < vr. We first 
handle the case that p\ lies in the spiral region of P2 ■ In this case we can connect the points by an 
inward going path from P2 to p\ using two spiral segments. By Equation [T] the length of this path 
is sec(a)(-R2 — Ri) < sec(a)D, which proves the claim for this case. 

Next we handle the case that p\ does not lie in the spiral region of p2- In this case we join 
the points using the right spiral through p\ and the left spiral through p2- Let p = (R,(j>) be 
the point where the two points first join. The summed length of the spiral segments is L = 
sec(a)(i?i + R2 — 2R), which we need to bound in terms of D 

We distinguish two cases. First assume the points have a distance of at most 3D/2 to the root. 
Then we obtain the connection between the points by simply connecting both to the root. Then 
L < sec(a)(R\ + R2) < 2sec(a)3D/2 = 3Dsec(a). Next assume R2 > 3D/2. From the discussion 
of Equation [2] above we know that L is maximized for i?i = R2. In this case we have eft = (02 + 0i) /2 
and therefore 

L = sec(a)(2R 2 - 2R 2 e 2tan ( Q >) = sec(a)2i? 2 (l - e 2tan ( Q > ) 

< sec(a)2.R2 777 — t\ = csc(a)R 2 {4>2 - <t>i)- 
2 tan (a J 

It remains to bound R2{4>2 — <fii) in terms of the Euclidean distance D of the two points. Observe 
that for given P2 (with R2 > 3D/2) and D the angle <p2 — 4>i is maximized if the line through the 
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origin and p\ is tangent to the circle of radius D around p2- Thus 2 — 0i is maximized if the angle 
formed by p%, p\ and the origin is tt/2. In this case 02 — 0i = arcsin(Z) / R2) ■ Thus, in general 

f D\ ( D \ D 

02 — 01 < arcsm — - = arctan — : < — : . 

\y/R%-D*J - 

Since R2 > 3D/2, we have \J R\ — D 2 > Rzy/l — 4/9. Plugging this into the above bound gives 
02 — 4>i 5; -D/ (i?2 -y/5/9). Now inserting this into the bound on L gives 

L < cac(a)D/y/E/9 < 3csc(a)£>. 
Combining the cases results in the claimed bound. □ 

Theorem 2 The length of the optimal spiral tree of a set of points is bounded by 

3[~log 2 re] max(sec(a), csc(a)) times the length of the minimum spanning tree of the set of points 

with the origin included. 

Proof. In the following we construct a spiral tree for which this bound holds. Let L be the length 
of the minimum spanning tree on the points including the origin. Let k = 3max(sec(a), csc(a)). 
Let C\ be a cycle through the points of length at most 2L (e.g., obtained by ordering the points 
based on a depth first search in the minimum spanning tree). We replace each edge of C\ by an 
inward going sequence of at most two spiral segments. This results in a cycle C[ of sequences of 
spiral segments of length at most 2kL. By taking either every even or every odd sequence we join 
pairs (possibly leaving the root unmatched) of nodes by spiral segments of total length at most 
kL. We repeat the construction on the join nodes (and possibly an unmatched point)) using the 
cycle C2 induced by the order given by C[. Again we form a cycle of spiral segments through the 
vertices of C2 ■ If we parameterize corresponding sequences in the cycles C[ and C' 2 by the angle 
from the origin, the sequences in cycle C 2 are closer to the origin than the corresponding sequences 
in C[ for all angles. Thus as consequence of Equation [2] the length of C' 2 is bounded by the length 
of C[ and therefore by 2kL. Thus as in the previous step we can join pairs of join nodes using 
spiral segments of total length at most kL. Next we construct C3 from C 2 in the same way and 
iterate the construction. After at most (log 2 n\ all nodes have been joined. The total length is 
then |~log 2 n\ kL as claimed. □ 



3.2 Relation with rectilinear Steiner arborescences. 

Both rectilinear Steiner arborescences and spiral trees contain directed paths, from the root to 
the terminals or vice versa. Every edge of a rectilinear Steiner arborescence is restricted to point 
right or up, which is similar to the angle restriction of flux and spiral trees. In fact, there exists a 
transformation from rectilinear Steiner arborescences into spiral trees. Consider the following trans- 
formation from the coordinates (x, y) of a rectilinear Steiner arborescence to the polar coordinates 
(R, 0) of a spiral tree. 

R =e x+y (3) 
=(y — x) tan(a) 

Assume we keep one of the coordinates x or y fixed. Using the spiral equation from Section [2] we 
see that the result is a spiral. More specifically, keeping x fixed results in left spirals and keeping 
y fixed results in right spirals. So that means that the above transformation transforms horizontal 
and vertical lines to right and left spirals, respectively (see Figure [9]). The transformation maps 
the root of the rectilinear Steiner arborescence to (1,0). Thus, to get a valid spiral tree, we still 
need to connect (1,0) to r. 
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Figure 9: A rectilinear Steiner arborescence transformed to a spiral tree. 

Lemma 4 The transformation in Equation transforms a rectilinear Steiner arborescence into a 
spiral tree. 

Unfortunately, the transformation has several shortcomings. First of all, the transformation is not 
a bijection, it is a surjection. That means we can invert the transformation, but only if we restrict 
the domain in the rectilinear space. But most importantly, the metric does not carry over the 
transformation. That means that it is not necessarily true that the minimum rectilinear Steiner 
arborescence transforms to the optimal spiral tree. Thus the relation between the concepts cannot 
be used directly and algorithms developed for rectilinear Steiner arborescences cannot be simply 
modified to compute spiral trees. However, the same basic ideas can often be used in both settings. 

4 Computing spiral trees 

In this section we describe algorithms to compute (approximations of) optimal spiral trees. First we 
show that it is NP-hard to compute optimal flux or spiral trees. Then we give an exact algorithm 
for computing optimal spiral trees in the special case that all spiral regions are empty, i.e. t{ ^ Tit- 
for all i 7^ j. Finally we give an approximation algorithm for computing optimal spiral trees in the 
general case. 

4.1 Computing optimal flux and spiral trees is NP-hard 

For the hardness proofs we will choose a = 7r/4. The reduction is from the rectilinear Steiner 
arborescence (RSA) problem [18] for spiral trees, and from the Euclidean Steiner arborescence 
(ESA) problem [T7] for flux trees. Shi and Su [18] proved by a reduction from planar 3SAT that 
the decision versions of the RSA problem and the ESA problem are NP-hard. Their reduction uses 
points on an 0(m x m) grid, where m bounds the size of the 3SAT instance. We can assume the 
grid to be an integer grid. Now if there is a satisfying assignment, the optimal RSA and ESA have 
an integer length K, while if there is no such assignment the optimal RSA and ESA have length 
at least K + 1. 

We can therefore state the problems for which they proved NP-hardness and from which we will 
reduce as follows. 

Instance: A set of integer points P = {pi, ■ ■ ■ ,Pn} m the first quadrant of the plane with 
coordinates bounded by 0(N 2 ); a positive integer K. 

Question (RSA): Is there a RSA of total length K or less? Otherwise the shortest RSA has 
length at least K + 1. 

Question (ESA): Is there a ESA of total length K or less? Otherwise the shortest RSA has 
length at least K + 1. 
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The basic idea is sketched in Figure [TUJ Assume we are given an . ^ _ H 

instance of the Euclidean Steiner arborescence problem with polyno- 
mially bounded coordinates. We translate the set of terminals by a 

large polynomial factor along the diagonal with slope 1 and place a i 
new root at the origin. If the bound on the coordinates is small (the | 
square in Figure fTU|) relative to the factor of the translation, then the 1 
angle formed by the line through any of the translated points and the 'i 
origin with the x-axis is "more or less 7r/4". A 7r/4-restricted flux 

tree thus behaves within this square "more or less" like an Euclidean Figure 10: Reduction. 
Steiner arborescence. For spiral trees we use the same setup but show 

that the distances on the spiral tree approximate the Li-norm. However, quantifying "more or 
less" precisely is technically rather involved and will be done in this section. 

To draw the connection from Steiner arborescences to flux and spiral trees we generalize the 
concept of RSAs and ESAs. For flux and spiral trees the angle a is bounded relative to the root 
while for RSAs and ESAs the angle that an edge can make with the x-axis is bounded (or with any 
given line through the origin). For ESAs the angle of an edge is in [0, vr/2], while for RSAs the angle 
is in {0, 7r/2}. We call a Steiner arborescence with angles in [/3, n/2— (3] a (> /3)-Steiner arborescence 
((> /3)-SA) and a Steiner arborescence with angles in {/3, 7r/2 — /3} a /3-Steiner arborescence (/3-SA). 
We do not restrict (3 to be positive but to — 7r/4 < j3 < 7r/4. In the following Steiner arborescences 
are typically not rooted at the origin. 

Now let -7r/4 < f3 < f3' < ir/4. Every (> /3')-SA is a (> /3)-SA but the converse does not hold. 
However, we can transform a (> /3)-SA to (> /3')-SA of similar length. Our transformation first 
transforms the whole tree and then connects the original points to their images under the first 
transformation. The transformation actually changes the location of the root. For the reduction 
we give this is not a problem because in the reduction we will have an additional root to which 
both the root of the original tree and the root of the transformed tree need to connect. 

Let rji = (— cos f3, — sin /3) and r\2 = (— sin/3, — cos /3). Let pi, . . . ,p n be points with pi = 1^771 + 
Vif]2-, (ui,Vi) € [0, B] 2 , where B may depend on n. Let T be a (> /3)-SA on pi, . . . ,p n with root 
Brn+Bri2. Let A = cos(/3+/3')/cos(2/3') and rj x = A(- cos/3', - sin/3') and r{ 2 = A(- sin/3', - cos ft). 
We transform T by the following transformation r : M? — > M 2 : Any point p = urj\ + vr}2 is mapped 
to q = urj^ + vrl 2 - We obtain a Steiner arborescence T' on p\, . . . ,p n with root Br][ + Br]' 2 by 
connecting pi to r{pi) by a line segment. Let A' = sin(/3' — 0)f cos(/3 + /?'). 



Lemma 5 V is a (> /3')-SA and \T'\ < \\T\ + 2\'Bn. 

Proof. If we ignore the connections between the p^s and r(pi)s the resulting transformed tree by 
construction fulfils the angle restriction and its length is A|T|. We therefore only need to show that 
the connections fulfil the angle restriction and that the length of any connection is bounded by A'. 
We have p% — r{pi) = Ui(i]i — !][) + v i(r\2 — f]^)- It therefore suffices to prove that rji — r][ and 7/2 — 1I2 
fulfil the angle restriction. Since 7/2 — r]' 2 is r/i — r][ mirrored at the diagonal with slope 1, it actually 
suffices to consider 771 — rj'i . 

Consider the triangle formed by the origin, 771 and r][ (see Figure [TT]) . (0, 0) 

We have 1 771 1 = 1 and |^| = A. By the law of sines A = sin 7/ sin 7'. 
This equation holds for 7 = tt/2 + (3 + /?', since then 7' = ir — ((3' — (3) — 
7 = 7r/2 — 2/3' and therefore sin 7/ sin 7' = sin(7r/2 + /3 + /3') / sin(-/r/2 — 
2/3') = cos(/3+/3')/cos(2/3') = A. On the other hand with 7 = vr/2+/3+ 
/3' we have is indeed reachable from 771 in a (> /3')-SA. The length of 
the connection is here the length of the third side of the triangle, which Figure 11: Lemma [5j 
is sin(/3' — /3)/sin7 = A'. More generally the length of a connection is 

bounded by 2B, that is B for each coordinate. Since we have n such connections the bound of the 
lemma holds. □ 
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In Lemma Owe have two summands, one depending on \T\ and one on n. Since the terminals lie 
on an integer grid and since every terminal has to connect to the tree, the length of the tree is at 
least of order n. 

Observation 1 If the terminals of a Steiner arborescence T have integer coordinates then n < 2\T\. 
Theorem 3 It is NP-hard to compute the optimal Eux tree of a point set. 

Proof. Given an ESA instance with root (0, 0) and with the coordinates x and y of any point 
(x,y) on the tree bounded by era 2 , we translate every terminal by 2n fc (l, 1) for a constant integer 
k > 2 specified later. We include the translated root in the point set but not as root. Instead 
we take (0, 0) again as root The shortest ESA is simply the originally shortest translated with one 
additional edge from (2n k , 2n k ) to (0, 0). Now, consider a point (c'n k + x, c'n k + y) with c' > 1 and 
< x,y < cn 2 . The angle of a line through this point and the origin with the diagonal of slope 1 
is bounded by /3 max = cn 2 /n k = c/n k ~ 2 . Now, restricted to such points every (^ /3 max )-SA is a 
flux tree with a = 7r/4, and every such flux tree a (> — /3 max )-SA. Also, every (> /3 max )-SA is a 
Euclidean Steiner arborescence, and every Euclidean Steiner arborescence a — /^max 

)-SA. Thus, 

if we show that (> /3 max )-SAs approximate (> — /3 max )-SAs well, this directly implies that flux trees 
approximate Euclidean Steiner arborescences well. More specifically, if we want to show that the 
length of the shortest flux tree approximates the shortest Euclidean Steiner arborescence up to a 
precision of 1 (so that we can make the distinction between K and K + 1) then it is sufficient to 
prove that a (> /3 m ax)-SA T' can approximate a (> — /3 max )-SA T up to this precision. 

By Lemma [5] and Observation Q] we get \T'\ < X\T\ +2\'cn 2 n < |T|(A + icX'n 2 ). Now, A = 
l/cos(2/3 max ) < 1/(1 - c/n k ~ 2 ) and A' = sin(2/3 max ) < 2/3 max = c/n k ~ 2 . Thus, 

A + AcX'n 2 < 1/(1 - c/n k - 2 ) + 8c 2 /n fc " 4 = 1 + o(l/n 4 ) 

for k > 8. Since \T\ = 0(n 4 ), this allows us to approximate the length up a o(l)-term. Note 
that we still need to connect the root of \T\ and the root of \T'\ to (0,0). The length of this 
connection is slightly different because the roots of the trees are different, but the difference is 
negligible compared to the difference of \T'\ and \T\. □ 

It remains to prove NP-hardness for spiral trees. 
Theorem 4 It is NP-hard to compute the optimal spiral tree of a point set. 

Proof. We use the same construction as above for spiral trees but we start with a rectilinear 
Steiner tree instance instead of a Euclidean Steiner tree instance. To adapt the reduction it suffices 
to show that within the relevant part of the tree, that is the part in [2n k ,2n k + cn 2 ], the length 
of a spiral segment between two points p, q is up to a small error the same as c"\\p — q\\\ for a 
suitable constant factor c" . The length of the spiral is D(p,q) = sec(7r/4)||»| — \q\\ = V2\\p\ — \q\\. 
Assume x q < x p and y q < y p . Let q' = (x p ,y q ). We have — \q\\ = (\p\ — \q'\) + (\q'\ — \q\). The 
difference \p\ — \q'\ measures how much the distance to the origin decreases while moving from p 
to q' . Let yo = y p — y q be the length of the line segment between p and q' and let a: [0,yo] — > M 2 
be this line segment parameterized uniformly. For a point u let ^{u) be the angle formed by 
the line through the origin and u with the x-axis. We have \p\ — \q'\ = Jq° cos ^y(a(u))du. Now 
7r/4-/3 max < 7(cr(u)) < 7r/4 + /3 max and therefore l/v^-^max < cos7(a(n)) < l/v^ + ^max- Thus 
IM — \q'\ ~ Vo/V2\ < c/n k ~ 2 . By the same argument we have that \\q'\ — \q\ — xq/\/2\ < c/n k ~ 2 , 
where xq = x p — x q . Therefore, 

\D(p,q) - \\p - qWx] = \V2((\p\ - \q'\) + (\q'\ - \q\)) - \\p - 

< \V2(x /V2 + y /V2 + 2c/n k - 2 ) - \\p - q\\i\ 
= 2v / 2c/ra fc " 2 . 
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Since the length of the RSA instance is in 0(n 4 ), the difference between measuring the length of a 
spiral segment versus taking the Li-distance of endpoints of segments is in o(l) for k > 4. Thus, 
computing the optimal spiral tree is NP-hard. □ 

To prove NP-hardness it was sufficient to consider one value of a, namely a = tt/4. Nonetheless, 
it is an interesting problem whether with results also holds for a given smaller a. We believe that 
the NP-hardness proof in [TH] can be adapted to /3-SAs and (> /3)-SAs for < (3 < tt/4. With this 
we could also generalize our result to smaller a. 

4.2 Optimal spiral trees with empty spiral regions 

Assume we are given an input instance such that ti ^ TZtj for all i ^ j. We give an exact polynomial 
time algorithm that computes optimal spiral trees for input instances with this property. 

Before we discuss the algorithm, we first give a structural result on optimal spiral trees for these 
special instances. Assume all terminals are ordered radially (on angle) in counterclockwise direction 
around r and are numbered as such. This means that the first terminal t\ is arbitrary and the 
remaining terminals t2, ■ ■ ■ , t n follow this order. First note that, for these instances, every terminal 
is a leaf in any spiral tree. That is because no terminal can be reached by another terminal, so no 
terminal can have incoming edges. More important is the following result. 

Lemma 6 If the spiral regions of all terminals are empty, then the leaf order of any planar spiral 
tree follows the radial order of the terminals. 

Proof. Assume this is not case, so that the leaf order skips leafs ti, . . . , tj, or in other words jumps 
from to ij+i- Pick any terminal tk with i < k < j. Let ir be the path in the spiral tree from 
tj_i to ij+i- Because the leaf order jumps from to ti+i, no leaf is connected to the outside 
(as seen from r) of ir. However, because tk ^ T^ti-u tk ^ ^-t j+1 and ti-i,tj+i ^ 7&t k , the path ir 
cuts through lZt k separating t/. from r. So the only way for tk to be connected to r is to cross tt. 
Contradiction. □ 

Corollary 1 If the spiral regions of all terminals are empty, then the leaf order of the optimal 
spiral tree follows the radial order of the terminals. 

Using the above lemma we can use a simple dynamic programming algorithm to compute the 
optimal spiral tree. We simply solve all subproblems that ask for the optimal spiral subtree for a 
sequence of terminals ti, . . . ,tj. We require that this subtree is contained in the unbounded wedge 
Wij from the radial line through ti to the radial line through tj (see Figure [12] left). Define pij as 
the intersection of and 5 f ~ (pu = U). As every internal node has exactly two incoming edges 
(Observation [2|), we split the subtree into two subtrees at every internal node. To compute the 
optimal spiral tree for a sequence of terminals ti, . . . ,tj, we simply compute the optimal way to 
split the subtree into two subtrees by trying all possibilities. We then connect both subtrees to 




Figure 12: Left: The wedge Wij for terminals U, . . . ,tj. Right: pij is the optimal join point. 
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Pij. Note that, by Lemma El we need to check only j — i ways to split this subtree. If F(i,j) is 
the length of the optimal spiral subtree for the terminals ij, . . . ,tj (contained in Wij), then we can 
perform dynamic programming using the following recursive relation. 



Note that we allow j < i, because we have a cyclical order. However, the value of k in the above 
equation must be between i and j in the cyclical order. The distance function D is defined as in 
Equation [TJ 

Lemma 7 The function F(i,j) describes the length of the optimal spiral subtree for the terminals 
U, . . . ,tj contained in Wij. 

Proof. We prove the lemma by induction. If i = j, then F(i,j) = is clearly correct. If i ^ j, 
then, by Lemma JUJ we compute the minimum of all possible splits for the corresponding subtree. 
Let this split be between and tk+i- By induction, F(i, k) and F(k + 1, j) describe the lengths of 
the optimal subtrees. We need to show that p^ is the optimal point to join the subtrees. For the 
sake of contradiction, assume the optimal join point is p'^. This point must be in the intersection 
of lZt t , TZtj and (see Figure [12] right). This means that p'^ £ TZ Pij - We can replace the edges 
Pik -» Pij and P(fc+i)j -» Pij by the edges p ik -> Pij, P(k+i)j -» Pij, and p^ -> p' iy Since p'^ must be 
closer to r than p^, it follows from the definition of D in Equation [1] that this operation shortens 
the tree. Contradiction. □ 

The length of the optimal spiral tree is not necessarily given by F(l, n), but it can also be any of 
the lengths F(i, i — 1) for 2 < i < n, so we need to compute the minimum of all these values. Note 
that there must be at least one wedge Wi^_i^ that contains the entire optimal spiral tree, so this 
will give the length of the optimal spiral tree. Using additional information we can also compute 
the optimal spiral tree itself in this way. From the definition of F(i,j), it is clear that the algorithm 
runs in 0(n 3 ) time. 

4.3 Approximation algorithm 

As shown in Section computing the optimal spiral tree is NP-hard in general. In this section we 
describe a simple algorithm that computes a 2-approximation of the optimal spiral tree. Note that, 
using Theorem [H this algorithm also directly computes a (2 sec(a))-approximation of the optimal 
flux tree. 

For rectilinear Steiner arborescences, Rao et al. |16] describe a simple 2-approximation algorithm. 
The transformation mentioned in Section 13.21 does not preserve length, so we cannot use this 
algorithm for spiral trees. However, below we show how to use the same global approach — sweep 
over the terminals from the outside in — to compute a 2-approximation for optimal spiral trees in 
0(n log n) time. 

The basic idea is to iteratively join two nodes, possibly using a Steiner node, until all terminals 
are connected in a single tree T, the greedy spiral tree. Initially, T is a forest. We say that a node 
(or terminal) is active if it does not have a parent in T. In every step, we join the two active nodes 
for which the join point is farthest from r. The join point p uv of two nodes u and v is the farthest 
point p from r such that p G 1Z U n 1Z V . This point is unique if u, v and r are not collinear. 




(4) 
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Lemma 8 The greedy spiral tree is planar. 




Proof. Assume there is a crossing in the greedy spiral tree between two u 2 q pv\ 

spiral segments, one between u\ and its parent v\, and another between 
U2 and its parent V2- Note that the intersection must be farther from r 
than both v\ and V2- But that means that the intersection must have been 
encountered while both u\ and v>2 were in W, so this intersection should be <* u i "-'^ 

a node in the greedy spiral tree. Contradiction. □ 

The algorithm sweeps a circle C, centered at r, 
inwards over all terminals. All active nodes that lie 
outside of C form the wavefront W (the black nodes 
in Figure [13]) . W is implemented balanced bi- 
nary search tree, where nodes are sorted according 
to the radial order around r. We join two active 
nodes u and v as soon as C passes over p uv . For 
any two nodes u, v G W it holds that u £ 1Z V . By 

t rout, j i j i Figure 13: The wavefront W. 

Lemma [8| the greedy spiral tree is planar, so we can ° 

apply Lemma [6] to the nodes in W. Hence, when C passes over p uv and both nodes u and v are 
still active, then u and v must be neighbors in W. We process the following events. 

Terminal. When C reaches a terminal t, we add t to W. We need to check whether there exists 
a neighbor v of f in W such that t G 7£„. If such a node u exists, then we remove v from W 
and connect u to t. Finally we compute new join point events for t and its neighbors in W. 

Join point. When C reaches a join point p uv (and u and u are still active), we connect u and v to 
Pmd- Next, we remove u and u from W and we add p uv to W Steiner node. Finally we 
compute new join point events for p uv and its neighbors in W. 

We store the events in a priority queue Q, ordered by decreasing distance to r. Initially Q contains 
all terminal events. Every join point event adds a node to T and every node generates at most 
two join point events, so the total number of events is 0{n). We can handle a single event in 
O(logra) time, so the total running time is 0(n log n). Next we prove that the greedy spiral tree is 
an approximation of the optimal spiral tree. 

Lemma 9 Let C be any circle centered at r and let T and T' be the optimal spiral tree and the 
greedy spiral tree, respectively. Then \C n T'\ < 2\C n T\ holds where \C n T'\ is the number of 
intersection points between C and T' . 

Proof. It is easy to see that \C Pi T'\ = \W\ when the sweeping circle is C. Let the nodes of W be 
ui, . . . , Uk, in radial order. Any node itj is either a terminal or it is the intersection of two spirals 
originating from two terminals, which we call uf and uf (see Figure fl~4"|) . We can assume the latter 
is always the case, as we can set u\ = m = if Ui is a terminal. Next, let the intersections of 
T with C be v\, . . . , Vh, in the same radial order as u±, . . . , ttfc. As T has the same terminals as T', 
every terminal uf and uf' must be able to reach a point vj. Let if and if be the reachable parts 
(intervals) of C for uf and uf, respectively (that is if = C fl 1Z u l and lf = CC\ H u r). Since any 
two neighboring nodes and iti+i have not been joined by the greedy algorithm, we know that 
if H lf_ 1 = 0. Now consider the collection Sj of intervals that contain Vj. We always treat if and 
if as different intervals, even if they coincide. The union of all Sj has cardinality 2k. If \Sj\ > 5, 
then its intervals cannot be consecutive (i.e. if, if, lf + i, lf+\, etc.), as this would mean it contains 
both if and ifi^ for some i. So say the intervals of Sj are not consecutive and Sj contains if and 
ih-li Du t not if (other cases are similar). T' is planar, so this is possible only if if C lf +l (see 
b igure US]). But then if and lf +1 are both in a collection Sji and we can remove lf +l from Sj, 
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Figure 14: Nodes ttj G W, terminals uf,u^ and intervals Figure 15: C 

I L I R 

while keeping the union of all collections the same. We repeat this process to construct reduced 
collections Sj such that the union of all collections remains the same and all intervals in a collection 
Sj are consecutive. As a result, |<Sj| < 4, and hence Ah > 2k or k < 2h. □ 



Theorem 5 The greedy spiral tree is a 2-approximation of the optimal spiral tree and can be 
computed in 0(n log n) time. 

Proof. The time bound is already mentioned above. For the approximation, recall that L(T) = 
sec(a) / °° \T n Cr\cIR, where T is any spiral tree and Cr is the circle of radius R centered at r. 
Using Lemma O we can directly conclude that the greedy spiral tree is a 2-approximation of the 
optimal spiral tree. □ 

The approximation factor is most likely not tight. Experiments for rectilinear Steiner arborescences 
show that the greedy algorithm often computes near-optimal arborescences [8]. 

5 Approximating spiral trees in the presence of obstacles 

In this section we extend the approximation algorithm of Section 14.31 to include obstacles. Given 
the similarities between spiral trees and rectilinear Steiner arborescences described in Section \3.2\ it 
makes sense to consider existing algorithms for rectilinear Steiner arborescences in the presence of 
obstacles. Unfortunately, the only known algorithm for this seems to have some issues. We discuss 
these issues in the next section. Then we give a new algorithm for computing rectilinear Steiner 
arborescences in the presence of obstacles. For a certain type of obstacles, this algorithm also 
computes a 2-approximation of the optimal rectilinear Steiner arborescence, although this does not 
hold for general obstacles. Finally we extend this algorithm to compute spiral trees in the presence 
of obstacles, again computing a 2-approximation for a certain type of obstacles. 

5.1 Ramnath's algorithm 

Ramnath [15] gives a 2-approximation algorithm for rectilinear Steiner arbores- 
cences with rectangular obstacles. He claims that the result extends to arbitrary 
rectilinear obstacles. But this is not the case. Consider the configuration of 
points and obstacles as seen on the right. The obstacles are L-shaped with the 
longer (vertical) side of the L being much longer than the shorter (horizontal) 
one. Between each consecutive pair of obstacles their is a terminal. What Ram- 
nath's algorithm does is to sweep a line of slope —1 starting at the root. During 
the sweep the arborescence is constructed greedily maintaining a minimal set of points (called cover 
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points) on the sweep line such that all remaining points can still be connected. Thus, in the begin- 
ning the algorithm has to decide whether to grow the arborescence to the right or upwards. From 
these two options the algorithm picks an arbitrary one, in particular it might grow to the right. 
But then on the arborescence will connect to each terminal by a connection corresponding to the 
longer side of the L-shape. By making the L-shape sufficiently long, the approximation factor for 
this configuration can be made worse than any constant, in particular two. 

Ramnath's paper also lacks the details to establish the claimed running time for rectangular 
obstacles. In particular the subdivision of a critical region (that is, a region that can be exclusively 
reached by one of the cover points) seems to assume that there is no obstacle strictly inside the 
critical region. However this case might occur and it does not seem straightforward to extend the 
algorithm to handle this case. Furthermore, the algorithm needs to compute the point at which 
the critical regions of neighboring cover points meet. This point is found by tracing paths from 
both of the cover points. The cost of this tracing step does not seem to be handled in the analysis 
and it is not clear how to account for it. 

5.2 Rectilinear Steiner arborescences 

We are now given a root r at the origin, terminals t±, . . . , t n in the upper-right quadrant, and also m 
polygonal obstacles B\, . . . , B m with total complexity M. We place a bounding square around all 
terminals and the root and consider the "free space" between the obstacles as a polygonal domain 
P with m holes and M + 4 vertices. We describe a greedy algorithm that computes a rectilinear 
Steiner arborescence T, the greedy arborescence, inside P. Our algorithm returns only a topological 
representation of T. This can easily be extended to the explicit arborescence, which, however, can 
have arbitrarily high complexity. 

As before we incrementally join nodes until we have a complete arborescence. This time we sweep 
a diagonal line L over P towards r and maintain a wavefront W with all active nodes that L has 
passed. If L reaches a join point p uv of nodes u, v G W, we connect u and v to p uv and add the new 
Steiner node to W. Our greedy arborescence is restricted to grow inside the polygonal domain P. 
If a point p G P cannot reach r with a monotone path in P, then p is not a suitable join point. To 
simplify matters we compute a new polygonal domain P' from P, such that for every p G P' , there 
is a monotone path from p to r in P. For now we simply assume that we are given P' and that it 
has O(M) vertices. 

To compute join points we keep track of the reachable region of every R(p,q) v 

node u G W, that is, we keep track of the part of L that can be reached 
from u via a monotone path in P' . As soon as two nodes u,v G W can reach 
the same point p on L, then p is the join point p uv and we can connect u 
and v to p uv . To compute the path between u and v and p uv , we need some 
additional information. Here our definitions follow Mitchell [13]. Given two 
points p,q G P' (with x q < x p and y q < y p ), let R(p,q) be the rectangle with p and q as corners. 
We say that q is immediately accessible from p if p and q are in the same connected component 
of R(p, q) n P' and this connected component does not contain any other vertices or nodes. The 
parent of a point p G P' is the rightmost vertex or node from which p is immediately accessible. 
The topological representation of the greedy arborescence stores only the parent information. 

The status of the sweep line L consists of three types of intervals: (i) free intervals: points that 
cannot be reached by any node in W, (ii) obstacle intervals: points not in P', and (iii) reachable 
intervals: points reachable by a node in W. The latter type of interval is tagged with the unique 
node in W that can reach this interval. We split the reachable intervals such that every interval 
has a unique parent. The intervals are stored by their endpoints in a balanced binary search tree. 
Initially, the status of L consists of one obstacle interval. We distinguish three types of events, 
which are processed in order using a priority queue. 

Terminal event. When we encounter a terminal ij, there are two cases. Either the terminal is in 
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a free interval or in a reachable interval tagged by a node u. In the latter case, we connect u to t% 
(using the parent information) and replace u by t$ in W. Also, we replace all intervals tagged with 
u by free intervals and merge them where possible. In both cases, we start a new interval for ij. For 
the endpoints of this interval, we trace the intersections between L and the horizontal and vertical 
line through ij. Note that we also split an interval, so we add three intervals in total and remove 
one. For every new interval (or merged interval), we add vanishing events to the event queue. 
Vertex event. When we encounter a vertex v, then v can be in any 
type of interval. If v is in a free interval, then we add an obstacle 
interval, where the endpoints of the interval trace the edges of P' con- 
nected to v. If v is in an obstacle interval, then we add a free interval, 
where the endpoints of the interval trace the edges of P' connected to 
v. Otherwise, v is in a reachable interval or at the endpoint between a 
reachable interval and an obstacle interval. In the first case, we need 
to insert an obstacle interval at v, as described above. In both cases 
we need to set the parent of v and insert a new reachable interval for v (with the correct tag). Also, 
we need to follow the edge or edges of P' connected to v. This can create free intervals. If one of 
the endpoints of the reachable interval of v directly moves out of P', we do not need to add this 
endpoint, but we can use the endpoint of the obstacle interval instead. Note that we add only a 
constant number of intervals. For the new intervals, we add vanishing events to the event queue. 
Vanishing Interval. If an interval I vanishes, then there are different cases depending on the 
types of the neighboring intervals I\ and 12- Note that / vanishes at a point p where two endpoints 
meet. If I\ and I2 are reachable intervals with different tags u\ and «2, then p is the join point for u\ 
and u%. We join u\ and 112 at p, as described in the terminal event. Otherwise, we need to remove 
one of the two endpoints. An endpoint of an interval always follows an edge of P' or a vertical or 
horizontal line through a node in W or a vertex of P'. If I\ and I2 are obstacle intervals or free 
intervals, then we can just remove both endpoints of /. If I\ and I2 are reachable intervals with the 
same tag, then we keep the endpoint that follows a horizontal line (this follows the definition of a 
parent given above). If I\ and I2 are of different types, then we keep the endpoint of the obstacle 
interval if one is present and otherwise we keep the endpoint of the reachable interval. Again, we 
add vanishing event points to the event queue for every interval for which an endpoint has changed. 

The algorithm terminates when L reaches r, at which point we have one node left in W. Using 
the parent information in the status, we connect the final node with r. To compute P' from P we 
simply run the sweep line algorithm in the opposite direction, tracing the "reachable region" of r. 
The points that border a reachable interval and either a free or obstacle interval trace out P'. 

Lemma 10 The greedy arborescence can be computed in 0{{n + M) log(n + M)) time. 

Proof. First we give a bound for the number of events. Clearly, the number of terminal and vertex 
events are bounded by 0{n + M). This also means that the total number of intervals is bounded 
by 0(n + M), as we add a constant number of intervals at only these events. At every vanishing 
interval event we remove an interval, so the total number of events is 0{n + M). Also note that 
every event can generate only a constant number of events. This also means that P' has complexity 
0(M), as we add vertices to P' only at events. It is easy to see that all events can be executed in 
0(log(n + M)) time, except when we need to change all intervals tagged by a certain node u to free 
intervals. We can do this in 0(n u ) time (by simple bookkeeping), where n u is the number of intervals 
tagged by u. An interval can only once be changed to a free interval. Merging two neighboring free 
intervals removes one interval, so we can charge these operations to the total number of intervals. 
Furthermore, the topological representation of the greedy arborescence contains only the relevant 
vertices and nodes to compute the paths between nodes. Every vertex or node can occur only once 
in this representation. So the algorithm runs in 0((n + M) log(n + M)) time. □ 
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If P has only positive monotone holes, then the greedy arborescence is a 2-approximation of 
the optimal rectilinear Steiner arborescence. A hole is positive monotone if its boundary contains 
two points p and q such that both paths on the boundary from p to q are monotone in both the 
x-direction and the y-direction. In the next section we prove this result for spiral trees. The same 
arguments can directly be applied to prove the same result for rectilinear Steiner arborescences. 

Theorem 6 The greedy arborescence can be computed in 0((n + M) log(n + M)) time. If P has 
only positive monotone holes, then the greedy arborescence is a 2-approximation of the optimal 
rectilinear Steiner arborescence. 




5.3 Spiral trees 

We now describe how to adapt our algorithm to spiral trees; we concentrate mainly on the necessary 
changes. We again compute only a topological representation of the output and refer to the spiral 
tree which we compute as the greedy spiral tree. The sweep line is replaced by a sweeping circle C. 
A simple balanced binary search tree is still sufficient to store the intervals, using special cases to 
deal with the circular topology. 

We need to replace horizontal and vertical lines by right and left r^SR(p,q) 
spirals. For a given node or vertex u, the endpoints of its interval on 
C follow the intersections of 5+ and S~ with C. Given two points 
p and q (q E 1Z p ), let SR(p,q) be the spiral rectangle between p 
and q. The spiral rectangle between p and q is bounded by the two 
paths (these are unique) consisting of exactly two spiral segments 
connecting p to q (this is exactly a rectangle transformed by the transformation in Section I3.2p . 
The point q is immediately accessible from p if p and q are in the same connected component of 
SR(p, q) Pi P' and this connected component does not contain any other vertices or nodes. 

There is one subtlety. If the left or right spiral of a vertex v directly 
moves out of P, we can ignore it, as before. However, at the exact 
moment that this is no longer the case, we do need to trace this 
spiral. For rectilinear Steiner arborescences, this can only happen 
at vertices. For spiral trees, this can also happen at most two spiral 
points in the middle of an edge e. A point p on e is a spiral point 
if the angle between the line from p to r and the line through e is 

I^i£fur6 16" S"dii*cl1 points 
exactly a. We hence subdivide every edge of P at the spiral points. v 

In addition we also subdivide e at the closest point to r on e to ensure that every edge of P has a 

single intersection with C . 

Neither the algorithm presented in Section \5. 21 nor its adaptation to spiral trees gives a constant 

factor approximation. But, if we restrict the types of obstacles, they give 2-approximations. For 

rectilinear Steiner arborescences we have to use positive monotone obstacles, for spiral trees spiral 

monotone obstacles. An obstacle is spiral monotone if its boundary contains two points p and q 

such that both paths on the boundary from p to q are angle-restricted. 

Lemma 11 Let P be a polygonal domain with spiral monotone holes. Then all points on a circle 
C reachable from a node u lie inside a single circular interval I u C C with the property that every 
point in I u n P is reachable from u. 

Proof. Consider a point p G I u D P. Let 7Ti and 7T2 be the paths from u to the endpoints of 
I u . Repeat the following until we hit either ir\ or tt2- Move along the left spiral through p going 
outwards (from r). When we hit a hole, simply follow the outline of the hole until we can follow 
the left spiral again. Because P has only spiral monotone holes, we eventually reach either 7Ti or 
7T2- Hence p must be reachable from u. □ 
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Theorem 7 The greedy spiral tree can be computed in 0((n + M) log(n + M)) time. If P has 
only spiral monotone holes, then the greedy spiral tree is a 2-approximation of the optimal spiral 
tree. 



Proof. Correctness and running time follow from Lemma [10] and the discussion in Section 15.31 
Assume that P has only spiral monotone holes and let T and T' be the optimal and greedy spiral 
tree, respectively. By Lemma [TT] we can represent the part of a circle C that is reachable by a 
terminal t as a single interval It- We can now follow the proof of Lemma [9] with these intervals to 
show that \C n T'\ < 2\C n T\. This directly implies that, if P has only spiral monotone holes, the 
greedy spiral tree is a 2-approximation. □ 
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